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7 I. INTRODUCTION 



8 From the mathematical point of view, many-body perturbation theory (MBPT) is repre- 

9 sented by a number of recurrence equations which permit to construct a total wave function 

10 of many-body system up to the fixed amendment. The higher order excitations are involved, 

11 the larger number of corrections is found. Due to interactions of the valence electrons, this 

12 number drastically increases for open-shell atoms. Over the past decades, it has become 

13 generally accepted^""- to build the wave operator which maps states in selected model space 

14 onto exact states of complex system Hamiltonian, upon the generalized Bloch equation, in- 

15 troduced by Lindgren et. al.-"-. Although, the practical application of this equation is far 

16 from being as mere as the theory is. The situation exceptionally changed, when symbolic 

17 programming tools appropriated sufficiently high level. Particularly, we exclude a compu- 

18 tational software program Mathematical . There exist several symbolic packages, written 

19 to evaluate the diagrams or algebraic expressions of MBPT-^--ii. Nevertheless, the specific 

20 features, necessary in computations, are known, as a rule, to the authors only. Therefore we 

21 developed the symbolic package NCopertors^, and the computational part of the present 

22 study relies on this package. The package has been tested by generating the terms of the 

23 first-order wave operator and the second-order effective Hamiltonian. 

24 The other two parts of this paper are devoted to the construction of model space and the 

25 reduction of generated terms. A good survey to various aspects of the problems that merge 

26 the selection of model space can be found in the reports of Lindgren, Kutzelnigg et.al.— 

27 On the one hand, we follow the traditional MBPT^ by excluding three types of the one- 

28 electron orbitals: core, virtual (or excited) and valence. On the other hand, we propose an 

29 algebraic investment that accounts for the operation of valence creation and annihilation 

30 operators on the model space in a strict manner (Sec. |Tl])- In the result, we formulate 

31 a precise statement which determines the amount of Fock space operators by the differing 

32 types of one-electron orbitals. That is, it ascertains the terms of the wave operator producing 

33 none zero contributions to the effective interaction operator. The theorem enables to simplify 

34 algebraic procedure, improving further on analysis (especially reduction) of terms. To apply 

35 the suggested formulation of PT, we display all obtained one-body and two-body terms of 

36 the third-order effective Hamiltonian in irreducible tensor form (Sec. I11I|) . The reduction is 

37 performed by using the angular momentum theoryi^ii^. The reduction scheme is presented 
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38 in a versatile disposition form, thus it is suited with the coupled-cluster (CC) approach, if 

39 some elementary replacements are initiated. Meanwhile, we show how suggested algebraic 

40 approach results to a compactness and easily accessible form of terms. The method is based 

41 on the construction of effective n-particle matrix elements, rather than the writing of each 

42 algebraic expression side by side its diagrammatic visualization, as, for example, in the 

43 classical works of Ho, Lindgren et. al.—'^^ as well as in the recently appeared papers^^i"— . 

44 II. THE SELECTION OF MODEL SPACE 

45 A. Hilbert space 

46 Suppose given a set X = {/ifc}^^ of orthonormal A^-electron eigenstates hk of the central- 

47 field Hamiltonian Hq with the eigenvalues Eq^. The configuration state functions (CSF) 

48 are characterized by the sets {U^ AkMk} of numbers: (i) the configuration parity 11^ = ±1; 

49 (ii) the irreducible representations of group G, where = LkSk if G = S0^(3) x SU"^(2) 

50 and Afc = Jfc if G = SU(2). The indices = M^^Ms^ or = Mj^ enumerate the basis 

51 for Afc. 

52 By^^, we let {,■,)-}{ '■ X x X — > M be a real function, called the scalar product on 

53 A^-particle Hilbert space Ti, where 



54 The infinite Hilbert space "H is assumed to be separable, thus there exist the linear forms 

55 G "H, the parameter e > and the integer such that — J^t—i Ckii)hk\\H < ^ for 

56 any M > and Cfc(z) G M, z = 1, 2, . . . , D. If particularly M = oo, then "^{i) will denote 

57 the eigenstates of the A^-electron atomic Hamiltonian H = Hq + V with the eigenvalues Ei, 

58 where V represents the perturbation. 

59 Along with the scalar product (or equivalently the norm ||, || = {,-,)), we also propose 

60 a unit operator on H. 

61 Proposition II. 1. The form 1 : — > H, expressed by 




n = {h 



{hk ■ hi)u = Ski,Wk,l G Z+}. 



(1) 




oo 




(2) 



k=l 



62 is a unit operator on N -particle Hilbert space H. 
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63 Proof. For the vectors hk, we immediately gather Ihk = Yl'i^i ^li^i ' hk)n = hk. In Ti, we 

64 may construct various hnear forms \& = Xlfeii '^fc^fc G "H, G M. Thus for the vectors we 

65 get = X^fcli Cfcl/ifc = ^- This proves the proposition. □ 

66 B. Orthogonal subspaces 

67 From the given set X of orthonormal functions (Sec. Ill Ap . we form a subset Y G X 

68 that is formed from the countable functions /i^ = 0^ VA; = 1, 2, . . . , d < oo. We insist the 

69 functions (pk to be identified by the sets of numbers {n^A^Mfc}. Here the configuration 

70 parity 11^ is identical for all 0a:, for all A; = 1, 2, . . . , c?. The representations are obtained 

71 by reducing the Kronecker products of irreducible representations A which in turn label the 

72 one-electron orbitals. These orbitals are represented by two types^: core (c) and valence (v). 

73 We also insist the subset Y to be complete by means of the allocation of valence orbitals in 

74 all possible ways. 

75 Since Y = {(f)k}i=i is the subset of vectors of Ti, it is sufficient to introduce a finite- 

76 dimensional subspace V (dim P = of "H by 

P = : ■ (l)i)n = Ski.yk, / = 1, 2, . . . , rf}. (3) 

77 Then the orthogonal complement = Q = 'HQV oiV is defined by 

Q = {eu = hd+k : {ek ■ ei)n = hi^k^ie z+}. (4) 

78 This immediately implies that 

{4>k-0i)n = \/k = l,2,...,d V/GZ+, (5) 

79 thus the orthonormal functions 6i form the complement Z = X\Y = {6i}'^i, and they are 

80 particularly characterized by the sets of numbers {Ilf A;M;}. Let us study the properties 

81 of 6i that are predetermined by the orthogonality in Eq. (jS]). First of all, Eq. is true, 

82 regardless of whether Uf = 11^ or not, as Z fl F = and E Ti (see Eq. ([T])). Secondly, 

83 Eq. does not exclude the functions 6i that could contain core and/or valence orbitals of 

84 0fc if nf 7^ n^. By Eq. ([1]), on the other hand, we are not confined to form the infinite set 

85 X of functions hk in any manner if reserving the formal conditions, presented in Sec. Ill A[ 
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86 Considerably, we insist the functions 6i to include the one-electron orbitals that are absent 

87 in all (pk, irrespectively whether Ilf = 11^ or not. These orbitals will be called excited (e) 

88 or virtual. 

89 C. The model functions 

90 From now on, we assume that the number D of the selected eigenstates \E'(i) of H equals 

91 to dim V = d. Then by Prop. III.lt 

d oo 
k=l 1=1 

d 

=$(i) + Q^(z), ^{i) = P^(z) = Cfc(0<^fc, (6) 

k=l 

92 where 

d oo 
k=l 1=1 

Cfe(z) = {4>k ■ = ■ Hi))H e (8) 

93 The functions will be called the model functions of V (see Ref.-). In Eq. ([6]), if replacing 

94 Q with Q = (fi — 1)P, then we simply get \l'(z) = where is called the wave operator. 

95 The procedure how to construct the generalized Bloch equation for f2 is well known and it 

96 will not be presented here. 

97 D. The improvement of generation of Hilbert space operators 

98 In this section, we wish to select the rules that allow to generate the terms of the wave 

99 operator as well as the effective interaction operator efficiently. As it will be demonstrated 

100 later, the rules under consideration significantly improve the computation of terms of higher- 

101 order PT. 

102 The wave operator is represented by the infinite series, where the first term is 1. The rest 

103 of terms Q^"^^ (m > 1) consist of the sum of n-body parts where n = 1, 2, . . . depends 

104 on m. We define fi„ (for any m, in general) by 
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Ol02---an/3l/32---/3n' 



(9) 




(10) 



105 In Eq. the summation is performed over all types (e, v, c) of the one-electron orbitals. 

106 Hereafter, we do not write in the sum the concrete values of orbitals (. . . Oj, . . . , . . .), only 

107 designating their type (a^ = Zj, (3j = Zj, z = e,v, c). The quantities 0L!ai...a„^i...i3,^ denote 

108 some structure coefficients, composed of the product of one-particle and/or two-particle 

109 matrix elements (energy denominators included). In accordance with Eq. ([7]), the Fock 
no space operators (creation) and (annihilation) are assigned as follows 



111 On the one hand, items (A)-(B) correspond to the definitions, presented in Ref.- (Eq. (13.3), 

112 p. 288). On the other hand, items (C)-(D) are strictly determined and they realize the 

113 arguments^ (Sec. 13.1.2, p. 288) that it is possible to create as well as annihilate valence 

114 electrons in V. It will be demonstrated later that these items are of special significance. 

115 Moreover, items (C)-(D) point to the definition of the complete model spaced by means of 

116 the allocation of valence states in 0^ in all possible ways (see Sec. IIIBp . 

117 The definition of f2„ conforms to the elucidation, enunciated in accordance with the 

118 generalized Bloch equation^ (Eq. (11.62), p. 247), as f2„ plays a role of operator, connecting 

119 Q and V spaces. (In diagrammatic representation i7„ includes open diagrams.) 

120 Here are the rules which establish the distributions of f3j providing none zero contri- 

121 butions of 0„(a/3). 



(A) a^P = 0, (C) OvP^O, 

(B) 4p = o, (d) 4p^o. 



(11) 



122 Proposition II.2 (^ p. 292). The operator P0n{af3)P ^ Va,/3 = v. 



123 Proposition II. 3. The operator Q0n{af3)Q ^ Wa,(3 = z. 



124 Proposition II. 4. The operator Q0n{a/3)P 7^ Va = e, v, V/3 = v, c. 



125 By Eq. (fTOj) . the self-adjoint operator 0^(a/3) = 0„(/3a). Thus, Proposition III. 41 may be 

126 reformulated in a distinct way. 
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127 Corollary II. 5. The operator POn{o:f3)Q 7^ Va = v, c, V/3 = e, v. 

128 The proof of Proposition III.2I is obvious, and it directly follows from Eq. ( fTTl) . 

129 Proof (Proposition III.3|) . It suffices to prove that Qa^Q 7^ 0, since in this case, the ex- 

130 pression QalQ = {Qa^QY 7^ 0. By Eq. (ITT]) and the expression P + Q = 1, we obtain: 

131 (1) QttcQ = Qac- But Qac = {alQY 7^ 0. (2) The operator Qa^Q = a^Q 7^ 0. (3) 

132 QayQ = [ttv, P] 7^ 0, where [■, ■] denotes a commutator. 

133 n 

134 The proof of Proposition III. 41 immediately follows from Propositions III.2tiII.3( and there is 

135 no point to present it here. 

136 Proposition |IL4] agrees with the statement of Lindgren^ (p. 292), if the identity a = /3 = v 

137 is neglected. The last condition, as a rule, is simply postulated to be false; otherwise, due to 

138 zero-valued energy denominators, the infinite terms f2„ are observed. However, we will study 

139 the operators (50„(vv)P on the Hilbert space "H in a more detail to show in which cases 

140 the rejection of the excluded condition a = /3 = v is true, since it has a direct connection 

141 to the properties of the model functions (see Eq. (E])) which form the subspace V (see 

142 Eq. OS])). 

143 Suppose for simplicity that n = 1. Also, let us mark (see Sec. IIIBI) 

cP, = /i(A^- A^- . . . Af-r,n^AfcMfc) , k = l,2,...d, 

144 where A^" is the uth electron shell in 0^; Tk denotes additional quantum numbers. Then 

145 by Eq. ([7]), the operator 



d 00 



QOi (vv) P = Y, ■ <P'k)H<Pl (12) 



k=l 1=1 



01 ^ a.,4,0. = (-l)E-\A^..+E-iiv..+^.+i^-,^^,^^ 

• • • ^fc j-iAf*'l'nfcjA^''''^^nA: j+iAf^ll' . . . nfc^A^^^IIfeAfcMfc) . (13) 

146 The parity of the iV-electron CSF 0'^ equals to 11^ = (— l)'''i+'^in^. It is assumed that 

147 and annihilate and create the zth and jth valence states of 0^, respectively. If = 11^, 



148 then, by the definition of V, obtained functions h{. . .) on the right hand side of Eq. ( IT3|) 

149 belong to the set F, and thus 

0'. = (-1)^^-^^^+^-^^^-^+'^^+^5a.,a,>.', A:' = 1, 2, . . . , 

150 Thus for i = j, we obtain a particular case k = k', if the shell A^'°' in 0^ is labeled by X^i = A^. 

151 In general, for some k, the functions (pk' are zeroes. Nevertheless, due to the completeness 

152 of the finite set Y (Sec. IllBI) . there will always be at least one function (pk with Xki = Ay. 

153 But {01 ■ (j)k')n = 0, VA;' = 1, 2, . . . , (see Eq. ([5])). This implies QOi (vv) P = 0. The 

154 generalization of studied example holds for any n. 

155 Corollary II. 6. The operator QOn{yy)P = 0, 'ifYll=i{^^i + ^fcvj 'is even for those k values 

156 which determine the electron shell UkYiX^^''^ in 0^ G V. 

157 (Obviously, the condition in Corollary II 1. 6 1 holds for POn{yY)Q too.) In other words, Corol- 

158 lary III. 61 represents an additional parity selection rule. The practical treatment of present 

159 rule in the study of triple excitations in CC approach can be found in Ref.^® (Sec. III-E, p. 

160 6). 

161 We can now summarize. The selected rf- dimensional subspace V of A^-electron separable 

162 Hilbert space "H is assumed to be constructed of the set Y of same parity configuration 

163 state functions (pk by allocating the valence electrons in all possible ways (complete model 

164 space). Additionally, in order to avoid the divergence of terms of the PT, we select the 

165 parity conservation rule in Corollary III. 61 to be true. The subspace V will be called the 

166 model space. 

167 Finally, let us define the effective interaction operator by^ (Eq. (15.5), p. 386) 

^ = = EE{^(^i + ^^)^nP}i, (14) 

168 where ^ = 1,2,3,4 denotes the number of contractions between the z-body parts (i = 1,2) 

169 of perturbation V and the n-body parts of Vt (for n + i — ^>Qi). On behalf of the definition 

170 of P, the associated operators Wn,^^ are generated by applying the following theorem. 

171 Theorem II. 7. // the n-body part of the wave operator is defined by Eq. then none 

172 zero terms of the effective interaction operator W on the model space V are generated by 

173 maximal eight types of the Fock space operators On{af3) for all n > 2. 
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174 Proof. We refer to Proposition 111.41 and Corollary III.6I None zero terms W^n,g are obtained, 

175 if at least for = 4, the operators {PsV2^nPs}4 7^ are derived. This means f2„ includes at 

176 least (n — 2) creation and {n — 2) annihilation operators, designating the valence orbitals. 

177 Possible allocations of and f3j (see Eq. flTU]) ) orbitals for all i, j = 1,2, . . . ,n are these: 



(1) 



/3, = V, 
i = 2,3 

I J = 



(2) 



n 



1,2, ...,n 



(4) 



«! = e, 02 = e , 



z = 3,4. 



(5) 



1,2, 



,n 

f 



/3i = c, (3j = Vj 
i = 1,2, . . . ,n 
J = 2,3,. . .,n 

Oti = Vi 

/3i = c, (32 = c' 

z = 1, 2, . . . , n 
j = 3, 4, . . . , n 



(3) 



(6) 



ai = e, ai= Vi 
h = c, (3j = Vj 
i = 2,3, . . . ,n 
J = 2,3,...,n 

ai = e, ai= Vj 
Pi = c, (32 = c', /3j 



= 2,3, 
I J = 3,4, 



n 



,n 



Ol2 



(7) 



^, = V, (8) <^ 



'3 — "J 

,n 



ai = Vi 
Pi = c, 
z = 3,4, , 
, j = 2,3,...,n 

179 Because of the anticommutation properties of creation and annihilation operators, the similar 

180 allocations of none valence orbitals hold for any selected i, j, not only i,j = 1,2. □ 



ai = e, a2 = e , cti 
Pi =c, /32 = c', f3j ■■ 
i = 3,4, . . . ,n 
, j = 3,4, . . . ,n 



V,- 



181 Theorem III. 71 also fits the effective Hamiltonian, given by the formula Hj^~^^^ 



'eff 



182 In this case, the structure coefficients to in Eq. are replaced with u^"^\ obtained from 

183 For m = 2, these coefficients will be displayed in the next section, where we examine 

184 a special case of the application of proposed formalism. 



185 III. THE THIRD-ORDER EFFECTIVE HAMILTONIAN 

186 In open-shell MBPT, the procedure to determine some fixed number i = 1,2, ... ,d of 

187 energy levels Ei of A^-electron atomic Hamiltonian H is addressed to the solution of eigen- 

188 value equations Heff^ii) = Ei^{i), where the model functions are determined in Eq. 

189 ([6]). However, in practical applications, the accuracy of effective Hamiltonian Heff is finite. 

190 In this section, we consider the third-order contribution to Hefj. 



191 The third-order approximation -f^gjj is represented by Eq. f|T^ . replacing fi„ with the 

192 second-order contribution fi*^^) = ^^=1 ■ Then 

2 4 min(2m,2n) 

^s^/=EE E = {^^"^^i^^^ie (15) 

m=l n=l ^=1 

193 The task under consideration is divided into two parts: (1) the determination of fin with 

194 n = 1, 2, 3, 4; (2) the construction of h^^^.c foT m + n — ^ = 1,2. 



195 A. The determination of terms of the second-order wave operator 

196 In the first part, the operators Qi? are generated in accordance with Ref.~ (Eq. (13.30), 

197 p. 302). The terms are computed by using the NCoperators package which is programmed 

198 upon Propositions III.2tiII.4l and Corollaries 111.51111.61 This part of computation is the most 

199 time consuming process. The generated terms are arranged by passing to Theorem III. 71 The 

200 coefficients cu^^j, cufj^p, u}^aKfiufi and w^^j^^^p.^-, located next to the operators Oi, O2, O3 and 

201 O4 in Eq. (Q, are treated as the effective one-, two-, three- and four-particle matrix elements 

202 on the basis of the product of accordingly same number one-electron eigenstates ip{na\a^a)- 

203 Here Aq, denotes the irreducible representation of G (Sec. Ill A|) . The basis index acquires the 



204 values nia = — Aq,, — Aq, + 1, . . . , Aq, — 1, Aq. Then the Wigner-Eckart theorem is applied to 

205 each matrix element, and the basis indices are excluded to the Clebsch-Gordan coefficients of 

206 SU(2). Despite of a large number of matrix elements uj^'^\ there are only a few of fundamental 

207 constructions to be examined; other elements are derived by varying the given ones. These 

208 S0(3)-invariant constructions are produced in Tabs. [Il lIIII All computed effective matrix 

209 elements u^'^\ necessary to form one-body and two-body terms of the effective Hamiltonian, 

210 are presented in an explicit form in Appendix |A1 The analysis of these elements is performed 

211 in the next part of computation. 

212 Let us study the structure J^^fi'^CP^iiaCis/i^ci^ ~ ^"a*) example. Here Eap...^ = £a + 

213 EjB + . . . + Et^, where Ea denotes the one-electron energy. The one-electron matrix element is 

214 given by V(;fi = {!f{n(;X(;m(;)-vJ^^!f{nfiXfimfi))sj = {n(^\c_m^\v'!^^\nfi,Xfimp) . The operator v acts 

215 on a single-particle Hilbert space S). Hereafter, the irreducible representation ri labels the 

216 one-electron operator v which befits to a second quantized form Vi in the generalized Bloch 

217 equation, written for VL^'^'^] the one-electron operator t>, located in ^2^^^ will be labeled by T2\ 
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TABLE L The multipliers for effective one-particle matrix elements of fi^^-* 



(ij^) Element: Expression 
(111) S^^inr^r) : 

X {TimiT2m2\Tm) 

(122) S^^in) : 

2(_1)A.+A, 5:^_(-l)Ac+A.^2M|) [,]V2 { n A, A, | 

(212) S'^i^{T2) : 

(223) 5,^- : 

4^A„A, [A.]-V2 E„ Ep,c(-i)'^-'''""i^ 

Xz{OXaX^XriXpUu)z{OXpXriXaX,^Uu) 



TABLE II. The multipliers for effective two-particle matrix elements of fi*-^-* 
(ijC) Element: Ex|")rossiou 
(110) D^p-p^iudr) : 

[r,,r2,u, d] /(-^^■^^mW^a.a.) ^^^^^^^^^ i^^^ 

r Aa A^ n 
X s Xjj, d 

Tl T2 T 

(121) Da^ijoiUuTi) : 

2(-l)A«-A,+A,+A.+n[[;]l/2 ^i^^(oA^ A^ApA^««) 

y / ^1 A« A<; "1 

^ \ A^ u U f 
(211) ^Up(^^^2) : 

R{S::i)Da0MUur2) 

(222) L>«^^p(nn) : 

4('_1)A;i+Ai7+«r^l-l/2 z(OXaXpXpX(Uu)z{QXpX(Xi?Xp,uu) 

(222) A^pi,i,{UU) : 

4(-l)^-^.[C/]V2 5^^^5^^^(_l)Ac+^JMi^ 
r At, A^ u ^ 

X < A,; d Xi? K z{OXaX(^XpXp,uu)z{OXpXfjXi^Xpdd) 

I. M Ap Ap J 
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TABLE III. The multipliers for effective three- and four-particle matrix elements of 17*-^-* 



Element: Expression 
(120) Tci3(;p,pfjiuTi) : 

R i'pl'^a) TaPCfivn (^1 '^2 ) 

(221) {-l)'''T^pciipr^{DdU) : 

4(_l)A,-A.+A,+t/[^]i/2 i_l)^[^]i/2,(0A,^/,ApA^nn) 

x.(oAAA,A,rfd){f^,'^^fJ{^,;J;;;} 

(220) Qal3Cpp.ufja{ud) : 

___4_(_-^^»+d+Ap+A.+A,+A.^(0^^^^^_^_^^ 

218 the operator v in Vi of Eq. (fT5|) will be labeled by tq. We assume that = (— l)'^('^i™-i)t) 

219 and = Vfi^. This implies 



/(nA^Ac) = e(riAcA^)/(riAcA^), (16) 
e(riAcA^) = (-l)^c-^^[ri]-i J2 (-l)'"-^(-i-^), (17) 

mi=— Ti 

/(riAcA,) = -^^[nc\c\\v^'\K\,], (18) 

220 where [x] = 2x + 1. The phase multiplier T is optional. Usually it is chosen to be equal to^ 

221 T(rimi) = Ti — mi. Then e(riA^A^) = (— l)'*'c^-'*M+n_ Particularly, if = rC^ represents 

222 the multipole momentum, then fiklf^li^) = f{kl(^lfi), since A^,^ = and the reduced 

223 matrix element [^(||C^||^^] 7^ 0, if + + is even. The state (p^UaXafna) denotes either 

224 2-spinor or 4-spinor, thus Eqs. (fT6l) - (fT8l) hold for both - none relativistic and relativistic - 

225 approaches. We assume that the one-electron Slater integrals are involved in the definition 

226 of /. The two-particle matrix element with tilde equals to Vp^^^^ = v^^^^ — Vp^^j^^, where 

227 the element f^^c/s ~ (^-^^i^^Ai^aAa'^alfi'ml^cAc^^C^/sA/!"^/!)- In general, the two-particle 

228 interaction operator gi2 acts on x i^. However, g'^ is reduced and it acts on irreducible 

229 tensor space S^^, obtained by reducing^^ (Sec. 2, Eq. (3)) S)^'^ x Sj'^'^. We also account 

230 for only scalar representations 7 = 0, and self-adjoint operators = 9i2- This implies 
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TABLE IV. The expansion coefficients for one-body terms of the third-order contribution to the 
effective Hamihonian 

^^^^^ (-l)^^-^^[ro]V2 ^^[A]i/2(^o"^oA M|AM) ((-1)^ Ee /(ToA.Ae)J]f +(A) { J a } 
-(-l)^Ec /(roAcA.)17ir (A) { - I A }) 

(212) 2(-l)^v-A. ^^^^^^ Ec(-l)''^'-'^ Eu z{OXA.X.Xa'Uu)n^l^{A)[u]y' { ^ J; } 
(-l)^N^/'EA,A,A(-l)^[Ai,A2,A]V2(^,^,AM|AAf)Ec(Ev'(-l)'^-'^7(ro^ 

xf)(?i+(AlA2A) A. A.A2 -Ee«(AeAvA2)/(ToAcAe)Cit(AlA2A) 
I TO A A J 





Av 




<^ Ac 


Av 


I}) 


I TO 


A 





2Ea,a,[Ai]1/2 Ecc' (a(AvA.A) Ev ?(0AeAc.A.,A,A2A2)oiJ+,(AiA2A) { ^1 ^ } 
(223) -"('^1'^^^) I(OAeA,, A,AeA2A2)J]fi:, (Ai A2A) { }) + 2 Ea,a. ("1)^^ i^]'^' 

X Ec { 1 1 1 } («(AvA.A) E^=v,e(-l)''''^'''"2^(0AvAcA,.A,,AiAi)17(^)+^,(AiA2A) 
-a(AiA2A)Eev'(-l)^^+^^'^(0AcAvAv'AeAiAi)17fj+(AiA2A)' 



2 EcC Em=v,c Ea, «(AcAc'A2) (^?(OAeAe. Av, A^,A2A2)Oi5+ ,,,,,(A2A2A0) 

(234) +Sa,A3a(-1)'^+^^+^[Ai,A3,A]V2(A3M3AM|AM)((-1)V' 

xI(0AeAe,A,,A.»A2A2)f)i^)+,,,,,(AiA2A3A) { i^;, t A^, } { t t xt } + "(^i^2A.) 
xg(0AeAe,A,»A,,A2A2)l^g)+,,,,,(AiA2A3A) { t A^ I { t I C }) 



''^p.ac^ = '^ap.jsc = '^pc.afi- Ref.^^, it was showed that f^a^/f may be constructed in two 

232 distinct ways: (i) reducing the Kronecker product (A^ x A(;) x (Aq x A^) (6-scheme); (ii) 

233 reducing the Kronecker product (A^ x A^) x (A^ x A^) (z-scheme). Then the two-particle 

234 reduced matrix element is formed in terms of either &(0A^AqA^A^7i7i) or z(OA^AQA^A^riri) 

235 coefficients^^ (Sec. 2, Eqs. (24), (30)) for 7 = 0. 

236 Usually the authors (see above cited works) better prefer 6-scheme, as it is more con- 

237 venient to produce the algebraic expression, written for a particular Goldstone diagram. 

238 Then a two-particle matrix element is denoted Vp,a.cp = Qfiacp Vf^^^cp = 9fia(;^j the cor- 

239 respondent reduced matrix elements are given by Xf^^(-p and Zf^^(^p (see, for example, Ref.- 
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240 (Appendix, Eqs. (Al), (A5))). Contrarily, in this paper we give priority exceptionally to 

241 the algebraic analysis of terms, and, along with Theorem IIL7t z-scheme is more preferable 

242 for the arrangement of their irreducible tensor form. Then w^q^^ (yfiaC^) "will be expressed 

243 by z{OXfiXaXpXcTiTi) (^OA^A^A^A^riri)), where 

^OA^A^A^-A^riri) = z{OXpX^X^X^T^T^) - a(A^A^ri)z(OA^A„A^A^riri), (19) 

244 where a(A^A^ri) = (— l)^c+^^+ri xhe studied effective one-particle matrix element is rep- 

245 resented by 



v-^^C^^ ^ (-l)^,«+™^-5,^(ri)(A„m„A^ -m^lnmi), (20) 

246 where S^^iTi) plays a role of the effective one-particle reduced matrix element. The quantity 

247 {XafnaXi^—mplTimi) denotes the Clebsch-Gordan coefficient of SU(2). Moreover, the form of 

248 Eq. fl20l) directly indicates that the one-particle operator, given by the product of Oi = aaO'^ 

249 and Eq. (!20!) . is simply equal to W^_^{XaX^)Sa^{Ti), where the irreducible tensor operator 

250 W^^{XaXi^) = [a^" X d^^Y^ is obtained by reducing Oi. The representation A^ designates the 

251 transposed annihilation operator dm^ = (— l)'*'^~"*^a^^^. The matrix elements of W^^ can 

252 be found in Kef.—'^. Here and elsewhere, it is considered if necessary that the summation is 

253 fulfilled over all given one-electron orbitals of marked type. However, only the sum running 

254 over the repetitive orbitals {(, ft in this case) will be written. 

255 In Tabs. [Il llllt the index i = 1,2 labels V^, while j = 1,2 labels ^f^- The operator 

256 R replaces orbitals in denominators. For example, the expression R{^^^a){^iyv — ^/3c)~^ 

257 reads {e^ — £a)~^- The quantities { ]^ ]^ ]^ } and < h h h ? denote 6j- and 9j-symbols. The 

258 elements which are found by making the three-pair contractions between V2 and Q'^^ vanish, 

259 if representations A^ 7^ A^ in W^i^XaX^). In this case, the orbitals (, p, r] are identical for 

260 all a, (3: ( = c, p = e, t] = v. Also, we mark off K^j^ = S^^jS'^^ by the summation 

261 parameter p = n^X^: (i) if yU = v, then K^^ = K^^, (ii) if = e, then K^^ = Kaj3i (iii) 

262 ii p = c, then we simply write K^^. The tildes designate that the direct and exchanged 

263 parts of a two-particle matrix element are involved. If given Da/3iiu{UuTi), then U marks v, 

264 represented by V2. Additionally, if given Dais/iuiUuTi), then u marks v, fitted to ^^2'^''. For 

265 both U and u, we write D^isfiuiUuTi). A similar argument holds for the rest of elements in 

266 Tabs. IITIIIIII These elements are also separated by the summation parameters. If ^ = 1, 
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TABLE V. The expansion coefficients for two-body terms of the third-order contribution to the 
effective Hamiltonian 

-N^/'EA,A[Ai,A]V2((-l)Aia(A.Avro)a(AiA2A)[Ai]V2(^o^oAM^ } 
(121) X Ee /(roA.Ae)ofj+ ,(AiA2A) { J; } + {-lr^>a{A^A,A)[A2]^/^{To - moA M|AM) 
X { Z t 1^ } /(roA.A,01^t?4(AiAiA) { - ^£ }) 



(211) 



(-1)V+a([A2]V2 ^^(-i)Ac2-(0A.A.,AeA.AiAi)Q(5+(A) { ^ ^ } " 
X ?(0AcA.A,,A.A2A2)f^S+(A) { ^ a^^, a^. }) 



a(AvAv'A2)[A2]-i/2 Ecc' I(0AcAe'A^,A^A2A2)17i5+,(AiA2A) + [Ai]-V2 

xEeE;.=v,e^(0AvA..A,.AeAiAi)O2+,,(AiA2A)+2(-l)^l+^^[Ai,A2]V2^^^^^J 

^^^^^ x[A2,^x]i/2 Eef(-1)^^+^^ Ee2^(0AvAcAeA^mn)J72+ (A1A2A) I „ A,, A A.j 
^ I Ac Av Ai_ A2_ J 

+ (_1)A.,+A. ^^„(_i)Ac+A,„~(0A,AcA,.A^,nn)o2i+ ^(AiA2A) | n a,, A a. ]) 



Ac Av Ai A2 



(-l)-o+Af[^^^^^^^]l/2^^^^^^^^ Ea2A3A^(-1)^[^2,A3,^]1/2(A3M3^^?|AM) 

x(romoA - M|^£.) ((-1)^-V+A[a,]1/2 Ea JAi]1/2/(toAcA,,0 



Av A2 A3 A Av Ai 1 

(132) x<^ A,, A,, Ul7^^!+ __,^(AiA2A3A) -(^/,va(Av'Av"A 

I A2 Ac A TO Ai X^i, ) ^ 
(2)+ M.T, A„T^^ _L ('_1^A,-,/+A„»+A2+A3+1?A f/'-T-^X \ ,Ao(2)H 



xJ^^Xwc(AlA2A3A)) + (-l)V+V,+A2+A3+^^^^j(^^^^^^,,)f^i2J^^^^^^^ 



f Ac 


A2 


A 




X <^ A,, 


Ai 






I Av 


A2 


A 


rn J ^ 



267 then the notations are similar to K^^ case. If ^ = 2 (see Tab. [TTl), then for Dai3p.u{uu), we 

268 write: (i) if C, p = c, then D = D; (ii) if (, p = v, then D = D; (iii) if ^, p = e, then D = D; 

269 (iv) if ^ = e and p = v, then D = D- For C = ^ and p = e, the similar triple doted 1) is 

270 considered. This is due to the symmetry properties of Vapjn). Finally, for Aa/3fiu{UU), we 

271 write: (i) if p = v, then A = A; (ii) if p = e, then A = A. In these cases, C = c. 
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272 B. The determination of terms of the third-order effective Hamiltonian 



273 By Proposition III.2t it follows that this part of computation requires significantly less 

274 time than the first one. Besides, the none zero terms of are derived in accordance with 

275 Theorem III. 71 which allows to reject a large amount of il^"^^ terms, attaching the zero- valued 

276 contributions. The operators h^^^.^ are considered by the formulas 

^!^i;€ = E<(^vA.)C^(A), (21) 

AM 

277 - for m + n — ^ = 1, and 



A1A2 

AM 

(3) (3)-\- 

278 - for m + n — ^ = 2. Each coefficient f)^„.^ is additionally expressed by the sum of f)^^.^ 

279 and f)^)jT^. The coefficients are presented in an explicit form in Tabs. HVirvTl The 

(3)— f3)+ 

2B0 coefficients f}.)„„.^ are derived from f)^^.^ by making the following alterations: 

2. (a) fi2+(A)^(-l)^^.+^-^+il]2-(A); 

- (b) ^^^(AiA^A) ^ (-l)^"''-+*^fiSr^(AiA2A); 

- (c) f^Sip,(AiA2A3A) ^ (-l)^"«'-'-+*^+^-^^+^fiSap,(AiA2A3A). 

284 Here Lap...c_ = + + . . . + A^. In addition, there holds one more rule: (d) each basis index 

285 (if such exists) in f)^)j^^, except for mo, is replaced by the opposite sign index. In (a)-(c), 

286 the indices M and M3 enumerate the basis for A and A3, respectively. The quantities QP''>^ 

287 are given in Appendix [Aj while 

f^Si(AiA2A) = ^^^^.(AiA^A) - a{K\,K{)^lf^%{K,K,K). (23) 

288 Particularly, the OS^'^^ are derived by replacing the one-particle and two-particle matrix 

289 elements v^p and Vaiifiv with and f^pa/j in order to obtain the standard form of the 

290 constructions that lay out in Tabs. HllIIII For example, the element Sap{riT2T) (Tab. [T]) 

291 is recognized from '^^Va^Vi^j^/ {e^ — e^) by excluding the S0(3)-invariant part. Thus, for 

292 a = V, /9 = c, H — e, we get Sy^ {tiT2t) which fits the definition of one-particle operator 
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TABLE VI. The expansion coefficients for two-body terms of the third-order contribution to the 
effective Hamiltonian (continued) 

[mnQ (A1A2A) 

2[A2]^/^Ec(-l)'^ EA,A.A3A(-l)^[A3,A]^/'((-imM3A-M|AM)[Ec'([Ai,Ai]V2 

x(_i)V+Ai ^^^^_^I(0AcAe'Av.A^.A2A2){Vf^S'tvc'c(AiA2A3A) - a(AvA^.Ai) 

xf^l?i'.c'c(AiA2A3A)} { I) 5 } I t V 1 1 - 5^^^^(-l)^c'+Ai 

'-'^'^^''lAaAAj 

X ?(0A,Ae.A,,A..A2A2)oJ)+,,,,,(AiA2A3A) { J } { M a_. a _ (_,)V+a, 

x[Ai,A2]l/2^^„ ^^__^(-l)^v"+A,,,+A,+A,~(o^^^^^_,,^^„^^^^)^(2H^^^ 

(233) X { 5, 5 J } 1 1' t t |] + (-l)^^+A^+^3+A[^^^Xi, A2] V2(A3M3A M|AM) EJnlV^ 



Av Ac A-/ Av A2 



X [(-1)^^ E^=v,e ?(OA.AeA^» A^.nn)0;;,^+ ^,,,,^(AiA2A3A) a,,, _ a, _ A, 

A^// A3 Ai A A^; A 

+(-1)"^' Ev" (Ec(-l)^^^(0AvAcAeAv,7X7x)f]i??Xw'c(Ai A2A3A) 

f M Av Ac A^/ Av A2 ^ 

Av" A2 Ai I - (-l)^i+«^-„(-l)V'l(OAvAcAv"Av"Uu) 



X 



[ Ae A3 Ai A A^/ A 

( u Av Ac A^/ Av A2 

xf^?/X"vv'c(AlA2A3A) <^ A,„ _ A2 _ Ai 

I A^„ A3 Ai A A^, A 



2<5a,a,(5ao(-1)^-+^-' Ecc'Aa (Ev"(Ev" a(Av"Av"Ai)2-(0AcAc'Av"Av"A2A2) 



x^^fvV'v''v'vc'c(AiAiA2A20) + EcA[A][a(AcAv''A2)[Ai,A2]"i/2~(0AeAc'A,"AeA2A2) 
x^^Svv'.'.c'c(A2AiAiA2A) - EA3AJA2, A3]V2(-i)A4a(A,A..Ai) { a. J } { Aa A. I 
xJ(0AcAc'A^,,AeA4A4){J^^JX,,^,^^,,(A2AiA3A4A) -a(Av'Av"A3)17^JJ,,,^,-,-^,^(A 
(244) + E^=v,c Ea[A][Ai, A2]-V2?(0AeAe.Ap.A^.A2A2)flg^„,,,,,,,,,(A2AiAiA2A) 

+ Ev"." Ea3A,a«(AiA3A.")[A][A2, A3]V2?(0AeA,,A,»A.»A4A4) { J J } { } 

X [^i"vv"v'v'vc'c 

(A2A1A3A4A) + o(AvAvA2)a(Av"Av"A3)0f-l^,^„-,-^,^(A2AiA3A4A) 
-a(A3A4Av){J^i'U,,„^,^,,,(A2AiA3A4A) + (-l)V'+A3a(A,Av,A2) 
xf^l?"v"v'.'.c'c(A2AiA3A4A)}] 



293 along with W^^X^Xc). On the other hand, another one-particle matrix element could be 

294 obtained from ^"m^m/s/I^m ~ ^a)- "^^^ ^^^^ element must be written in a standard form 

295 (—1) "^^v^i^Vf^a/ i^a — ^fi) to arrange the element Sp^{tiT2t) correctly. For concrete values 

296 a = V and /9 = c, the orbital /x equals to c only (see Proposition III.4p . and the element is 
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297 denoted Scv{tiT2t). In a tensor formalism, we gain the opposite sign of the basis index m 

298 in W^{Xj.c)- 

299 The diagrammatic interpretation of h^^~^ can be clarified as follows. The diagrams of 

300 the third-order effective Hamiltonian -f^e//' assembled in /i^^.^, are derived by contracting 

301 the perturbation Vm with iln^^ , where fin''" includes: (i) the folded diagrams; (ii) some 

302 diagrams, obtained by contracting the core orbitals in Wick's series; (iii) the diagrams, 

303 acceded to the reflection of a number of diagrams in fin about a horizontal axis. Note, 

304 ^apcpfii^w ^ ^ojcp/ii^'?^' ^^^^^ obtained from the disconnected diagrams V2^^2^ (see Tab. 

305 HID. 



306 C. The estimation of terms 

307 The evaluation of the number of computed one- body terms (Tab. Mi) of H^ff is presented 

308 in Tab. IVIIt where denotes the number of direct terms in h^^^^, while denotes the 

309 number of direct terms in h]^^.^, if the one-body interactions v'^^ {i = 0,1,2) are absent. 

310 Totally, there are computed 188 + 70 = 258 direct one-body terms of and 72 + 20 = 92 

311 direct one-body terms of h^^^.^, including the two-particle interactions (7° only. For instance, 

312 Blundell et. al.^^ calculated 84 diagrams contributing to the third-order mono-valent removal 

313 energy. In our considerations, their studied energies: a) E^^^-EP; b) and c) 
3u (see Ref.il (Sec. II, Eq. (8))) denote the matrix elements of terms in h^^,, h^l 

315 and h^2i-2i if 9^ represents the Coulomb interaction. 

316 The estimation of the amount of two-body terms (Tabs. IVTlVIl) that contribute to H^^^ 



317 is provided in Tab. IVIIII There are 217 + 82 = 299 direct two-body terms in h^^^c and 



''mn;^ 

318 125 + 4:2 = 167 direct two-body terms including the two-particle interactions only. In 

319 their study of beryllium and magnesium isoelectronic sequences. Ho et. al.- calculated 218 

320 two-body diagrams of the third-order perturbation. Analogous disposition to account for 

321 the two-particle interactions only, can be found and in other works^i^>2^. Additionally, most 

322 of them do not account for the folded diagrams. 

323 Meanwhile, the expressions in Tabs. HVllVIt obtained by exploiting the properties of 

324 proposed model space V (Sec. [Tll), have their own benefits: 

325 1. The third-order contributions to the effective Hamiltonian H^ff are written 

326 in an operator form providing an opportunity to construct their matrix elements 
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TABLE VII. The amount of one-body terms in MBPT(3) 







d+ 




d~ 




1 3 


n 


O 


n 


(122) 


37 





18 





(212) 


14 


2 


2 





(223) 


67 


34 


29 


2 


(234) 


57 


36 


18 


18 


Total: 


188 


72 


70 


20 



327 efficiently. Namely, the irreducible tensor operators, labeled by the representa- 

328 tions A (see Eqs. (l2Tl)-(l22l)). are written apart from the projection- independent 

329 parts. These angular coefficients include the structure coefficients fi*^^^^, mul- 

f jl 32 js J4 1 

330 tiplied by the 3nj-symbols. Particularly, the coefficients < h h h h \ 

I fci k2 kz ki J 

f ki k[ k k' k2 k'2 ~> 

331 and < Pi p P2 v (see Tabs. IVllVip . denote the 12j-symbol of the ffist 

l Jl j'l j j' h j'2 ' 

332 kindi^ (Sec. 4-33, Eq. (33.17), p. 207) and the 15j-symbol of the third kind^^ 

333 (Sec. 4-20, Eq. (20.3), p. 112). 

334 2. The form of the expressions in Tabs. llVtiVIt allows to evaluate the contri- 

335 butions of n-particle effects in CC approach. This is done by simply replacing 

336 fi*^^)^ with i7„ {n = 1,2,3,4). By generally accepted labeling, such replacement 

337 leads to the transformation gap,„(^ — )■ Pai3...c,i where p denotes the valence singles, 

338 doubles, triples, quadruples amplitude. 

339 3. Obtained terms of the third-order perturbation include, in addition, the one- 

340 particle operators v''^ {ji = 0, 1, 2, . . .) that represent the magnetic, hyperfine, etc. 

341 interactions. Moreover, the general expressions also ffi none relativistic as well 

342 as the relativistic approaches. These effects are embodied in ^(OAaA/^ApA^riri) 

343 coefficients. 

344 4. The reduction scheme allows to write a large number of terms in a concise 

345 form. This feature becomes evident especially clearly when in contrast the terms 

346 are written side by side their diagrammatic representation. 

(3) 

347 Maintaining the completeness of the present discuss, we note that -f^e// ^^^o includes the 
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TABLE VIII. The amount of two-body terms in MBPT(3) 





d+ 


d+ 


d- 


d- 


(121) 


20 





10 





(211) 


13 


2 


3 





(222) 


64 


32 


31 


4 


(132) 


20 


16 


10 


10 


(233) 


75 


50 


28 


28 


(244) 


25 


25 






Total: 


217 


125 


82 


42 



348 terms h^^i-e, '^ith m + n — ^ = 0,3,4,5. For example, the coefficient [)22.i (-E'iAiii^2A2A) reads 

[)g+(EiAiE2A2A) = (-l)'^"+'^'+'^"+^^+^[Ai,A2]'/'5](5]a(Av.A,A2) 

Ai A2 

X [Ei,E2,Ai]i/2^2^(0AeAvAv'AvMM)l^t'l+,^„(AiA2A) 

cu 

f Av A2 A2 1 f A2 Av A2 1 f A,, -El 1 r \-„ \-, E2 \ 
^ \ \.„ « Ac J \ Ai A Ai / \ Kn Ai Ai J \ Av A2 u j 

+ (_i)i^i+E2^^(0AvAv,AeAv£;i£;i)n2+_„(Ai£;2A) 

e 

f Ai Ai A^ 1 r Ai A2 A n ^241 

^ \ Ae £1 A,» / \ B2 Ai Av ' y^^^ 

349 and it forms the three-body operator h'^2-i along with the irreducible tensor operator 

[[ly^^AvAvO X a^'l^i X [W^^^(Av"Av') x d^'f']^. (25) 

350 However, the triple and higher-order effects are not covered by the examination of this paper. 

351 IV. CONCLUSIONS 

352 We present an algebraic technique to evaluate the terms of MBPT. The method relies on 

353 two main circumstances: the strictly determined operations of the Fock space operators on 

354 the vectors of the orthogonal subspaces of given separable Hilbert space, and the specific 

355 reduction scheme of terms of the PT. The aspiration to determine the behavior of the 
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356 creation and annihilation operators is motivated by the fact that in higher-order perturbation 

357 theories, a huge number of terms is generated and, particularly, not all computed terms of the 

358 wave operator attach none zero contributions to the terms of effective interaction operator. 

359 Therefore the rules that allow to predetermine these valuable terms become meaningful. 

360 Meanwhile, in the PT, another no less important procedure is to work up the generated 

361 terms in order to calculate their matrix elements efficiently. If going on the traditional 

362 route, when each term is expressed side by side its diagrammatic representation, then once 

363 again one deals with the redundancy of terms and the procedure of reduction becomes 

364 troublesome. In this paper, however, we suggest the reduction scheme which is of versatile 

365 disposition. That is, the generated terms are combined in groups related to the different 

366 types of one-electron orbitals. Afterwards, we construct the effective matrix elements and 

367 apply the Wigner-Eckart theorem. In the result, we get the irreducible tensor operators and 

368 the projection-independent parts that are located apart from each other. This means the 

369 form of reduced terms becomes universal, only the S0(3)-invariant parts that embody the 

370 dynamics of studied physical interactions are changed. Since the method of reduction is 

371 applied to the third-order MBPT, we obtain maximum four-body effective matrix elements 

372 that are derived from the four-body parts of the second-order wave operator. 

373 Finally, it is worth to mention that obtained symbolic preparation of terms of the 

374 MBPT(3) can be extremely simphfied if apphed to some special cases of interest. Mathe- 

375 matically, this means that some of the irreducible representations drawn in the 3nj-symbols, 

376 in most cases would be simply equal to zero. Particularly, the 12j- or 15j-symbols would 

377 become the ordinary and widely used 6j- or 9 j -symbols. Moreover, if the one-electron inter- 

378 actions, characterized by the none scalar representations, are not taken into account, then 

379 the irreducible one-body and two-body operators of the effective Hamiltonian are simply 

380 scalar operators. Thus their matrix elements become even simpler. 
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381 Appendix A: SO(3)-invariant peirt of the second-order wave operator 



382 One-body part. 

J7g+(A)(£e-£;.) 

~ ~ (Ala) 

= SAr[Si^c{TlT2T) + S ;j,c{tiT2t)] + Sati [Sficin) + Sf^dn)] 
+5At-2[5'^c(t2) + S'fj,ciT2)] + ^AoSf^c, 

42)- (A)(£e - e^) = SArSc^{riT2T). (Alb) 

~ ~ (A2a) 

= SATSev{TlT2T) + ^Ati [^ev (ti ) + Sevin)] 
+SAT2[S'ev{T2) + <S''ev(T2)] + (^AO'S'ev, 

ni^-{A){s^ - ee) = 5Ar[^ve(riT2T) + ^ve(Tir2T)]. (A2b) 

383 Two-body part. 



(A3a) 



^S'tc'(AiA2A)(£ec'-£^^') 

= SaiuSa^ciSatD^^/cc' (udr) + (5AiC/(JA2n'5ATi [^/i/i'cc'(^^''"l) 
+ |^W'cc'(f^«''"l)] + ^SAiuSA2uSAT2 2^tJ.'ficc'iUuT2) 
^D' ^i^,{UuT2) - ^AiA2<^Ai«^Ao[i-D^/i'cc'(^*'") 

{uu) + Ann'cc' {uu) + ^nn'cc'{uu)], 

^^2'cc'(AiA2A)(£cc'-£mm') 



= -5'^AiC7<5A2«'^Ari-2c'w'(wi7ri)5c'c;^^'(u?7ri) (A3b) 

+ |'^Ai!7<^A2«'^AT2^c'cAtAt' (^if^ T2) [^'c'cmm' (^^ '^2) 
+-D'c'cM/.'('"^'^2)] + 35AiA2'^Ai«(^A0-2cc>/.'(")-^cc'mm'('"^)- 
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J^f)+(AlA2A)(£cc'-£ev) 



+ h^vciJ,^'{uUTi)} + Zcyfj,'i^{uUTl)Dcyfj,',^{uUTl)] 

(uC/r2) 

{uUt2)} {uUt2)D' {uUr2)]. 



(A4a) 



= -5Ai„5A2d'^Ar^vccc'(^iC?''")I?vccc'(^iC^''") + (5AiC/5a2U<^Ati 
X [-Devcc' (C^'UTi ) + Zvec'c (f^ UTi)Dyec'c{U UTi ) 
+ 5-Devcc'(t^^tn) - ^Zyecc'iUuTi)Dyecc'{UuTi)] 
+ kSAiuSA2uSAT2[-D' evcc' iUuT2) + Zyecc'{UuT2) 
xD'yecc'{UuT2)] - ()^AlA2'^Al«<^A0[Acvcc'(^i^i) + Acvcc' 
-2'vecc' (■u){Avecc'(™) + Avecc'(^^)} + ^evc' c{u) D eve' duu)] , 



^evcc'(^1^2A)(ecc' -£ev) 



= -SAiuSA2U^ATiZc'ccviuUTl)Dc'ccviuUTl) (A4b) 

+(^Ai C/'^A2U<^Ar2 -^c'cev {uUt2) [-D'c'cev {uUt2) 

+D'^f^^y{uUT2)\. 

^S4(AiA2A)(ec.-eMM') 

= &Alu5A2d8^.T^^le[Dee'cv{udT) - Zee'vc{udT)Dee'vc{udT)] 

+SaiuSa2uSati [D^h'cv {UuTi ) - {UuTi ) ( A5a) 

X-Djit/i'vc(t^^^n)] + 5aiu5a2uSat2Z^i p,cy{UuT2)D' ^1 ^^{UUT2) 

+5A^A2^^lu5^o[Z^_i^_i'wc{u){/^^J.^J,'wc{uu) + /^^f^iyciuu)} 

^^2';.(AiA2A)(ec.-e^^o 



(A5b) 
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J^git(AlA2A)(£cv-£ev) 

= ^Aiu^A'zd^ATiDcvvciudT) + Z^ccv{udT)D^ccv{udT)] 
+^AiU^A2uSAti [Dcvvc{UuTi) + Dcvvc{UuTi) - Z^evciUuTi) 

X{Dy^yc{UUTl) + Dyey^iUUTl)}] + 6 5 A^J At^ i^^^) 

X [-D'eYYc{UuT2) + Zyeyc{UuT2)D'yewc{UuT2)] 
-5AiA2'^AiM<5Ao[Aevvc(^i^^) + Aevvc('"'") + DeYwc{uu) 

-Zyevc{u){Aycvc{uu) + ^vcvc{uu)} " ^cvcv (^i) { Agvcv (wu) 
"l"Aevcv 

(uu)} + Zyecv{u){^yecv{uu) + Ayecv (■"■«)}], 



i^£W'(AiA2A)(evv'-^ev) 



= - SaiU^AzuSati -Zy'vev (uUti) [-Dy' vev ('"t^Ti ) 
+Dy>yev{uUTi)] + 6A^u^A2uSAT2^v'vev{uUT2) 
xDy'yev{uUT2) - 5^AiA2<5Aiu<5A0-2^vv'ev(^^) 
xi)vv/ev(wS). 



(A6b) 



fi2k(AiA2A)(ecv -Se 



= SA^uSA2uSATi[2ycev{uUTi){Dycey{uUTi) - Dycey{uUTi)} 
-ZcyeviuUTi)Dcyev{uUTi)] + (5Ai!7(^A2W<5Ar2 [--^vcev (^^C^T2) 
XD'ycey{uUT2) + Zcyev{uUT2){D' cvev{uUT2) 
-b'^y^y{uUT2)]\. 

J^2i"t'(AiA2A)(e,,.-eev) 

= '^AiC/<^A2u'^Ari [-Oevvv' (^^ UTi) + 2^vevv' (^^ '"''"l) 

x{by^yy>{UuTi) " i yev v' ( C^-^Tl ) } ] + | <^Ai C/f^Aau (^At^ (A7a) 
X [Zvevv' (t^ Ur2)-D'vevv' (t^ ^^7-2) - -D'evvv' (C^'Wr2) 

-b'eyyy'{UuT2)] + (5ai Aa '^Ai«<^AO ^cvw' (^i^^) 
-Aevvv' + Zyeyyi{u){Ky^yyi{uu) + Avevv' (if") }] , 



(A7b) 
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Q2)vv'(^1^2A)(evv'-eee') 

{UuTi)] + |5Ail/'^A2M^AT2^e'ew'(^'"''"2) 
xD'c'cvv'{UuT2) - 5AiA2'^Aiu<Jao[|{^-Dcc'vv'(^^) 
-Zee'v'v(w)i>ee'v'v(Sli)} + Aee'vv'(^^) + ^ee'vv' i 



nfX-,{AiA2A){e^^' - See') 



(A8a) 



^ (A8b) 

X-C''v'vee'(Mt/r2) + j()^AiA2<JAiM'^A0^vv'ee' (^i) 
X [Dvv'ee' - -Dvv'oc' (uu)] ■ 

384 In Eq. (|A3ap . the quantities 3£^^'cc', 2)^^'cc', differ for distinct one-electron orbitals /x = v, e. 

385 If ^ = V, then Xvv'cc' = D^v'cc', Vw'cd = D^v'cc' and 3vv'cc' = -^vv'c'c(^i)-Dvv'c'c; if M = e, then 

386 Xee'cc' = -^ee'cc' 5 '^ee'cc' = -^ee'cc'5 3ec'cc' = -^ee'cc' • 

387 Three-body part. 



^wt"vc'c(^1^2A3A)(evc'c - evv'v") 

= |'^AA3'5MAf3[?'vv'v"vc'c(AlA2A) + |rvv'v"vc'c(AlA2A) 

- Ea,A Jvcvc'(A A2AA2Ai){Tvv'v"c'vc(AiA2A) (A9a) 

+Tvv'v"c'vc(AiA2A)}] + {-l)^''+^c'Y,,^^,{AiA2AA3Ti) 
x7Vv'v"c'cv(Airi) + i Xlu^lVAv-u) 

X /vv'v"vc'c 

(AiA2A3riAu)T' (^^7-2), 



^w'v"vc'c(^1^2A3A)(evc'c - ^vv'v" 



= 5<5AA3'5MAf3fl(-^v'Av"Ai)[a(AcAc'A2)Tvc'cvv'v"(A2AiA) 

+ EA2Aa(AcAc'A2)Jvc'vc(AA2AA2Ai)rcc' vvv V 

(A2A1A)]. 



(A9b) 
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^vv'evc'c(^1^2A3A)(£vc'c - £vv'e) 

= (^AAg <^MM3 [ 5 { ~ <^ (''^c Ac' A2 ) Tvy/ evcc' ( Al A2 A) 

+ 5T'vv'evc'c(AlA2A)} - a(AeAv'Ai){Tvev'vc'c(AlA2A) 

+ |rvev'vc'c(AlA2A)} - ^^^^ Jvcvc'(A A2AA2A1) 
x{i(rvv'ec'vc(AlA2A) + f vv'ec'vcCAlAsA)) 

+a(AeAv'Ai)a(AcAvA2)(rvev'c'cv(AiA2A) (AlOa) 
+rvev'c'cv(AiA2A))}] + (-l)^-+Vy^^-^,(AiA2AA3ri) 

X Tyv'ec'cv 

(Ain) + ia(AeAvA2)F;^,,^(riAiAA3A2) 

xrevv'vc'c(A2ri) + E„[a(AcAv'A2) 

X/vv'evc'c(AlA2A3riAu)rcvv'c'cv(^iTl) + ^a(Av'Avn) 
Xi'vv'evc'c(AlA2A3T2Au)f'evv'cvc'(^*^2) + (-1)^-"^^ 
Xa(AeAv'Ai)7vev'vc'c(AlA2A3T2An)r'v/ 

evcvc' 



^l?'evc'c(^1^2A3A)(evc'c " 



= (JAA3<^MM3a(AeAv'Ai)[a(AcAc'A2)rvc/cvv'e(A2AiA) 

_ _ _ (AlOb) 

+ 1 Ea2a{<^vcvc'(A A2AA2Ai)rc, 

cvvv'e 

(A2A1A) 

+ (_l)Av-Ae /;^,^_^,^(AiA2AiA2AA)(re',w(^2AlA) 

-a(AcAvA2)rc/ 

vcevv' 

(A2A1A))}]. 
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^ee'vvc'c(^1^2A3A)(£vc'c - ^ee'v) 

= 5aA3<5mM3[H^cc'vvc'c(AiA2A) + |Tcc'vvc'c(AlA2A)} 
+a(AvAe'Ai){a(AcAc'A2)reve'vcc'(AlA2A) 
-5^eve'vc'c(AlA2A)} + ^^^^^ Jgcvc' (A A2AA2A1) 
x{i(a(AcAvA2)Tee'vc'cv(AlA2A) - ree'vc'vc(AlA2A)) 

(Alia) 

+a(AvAe'Ai)(rcvc'c'vc(AlA2A) + rcvc'c'vc(AlA2A))}] 
+<5AiA25A3n'5A0ree'vvc'c(Airi) + (-1)^-=+^ 
xyecvc'(AlA2AA3ri)Tee' 

vc'cv 

(Ain) + E„[a(Ae'Ac'A2) 

^ -^ee'vvcc' (AlA2A3riAu)rvee'c'cv(^^T-l) + a(Ae'Av'u) 
2 -^ee'vvc'c 

(AiA2A3T2Au)T' 

vee'cvc' 

+ ( _ 1) Ae- Av 7g^g, ( A2 Asrs An) f 'e' vecvc' (?^r2 ) }] , 
^2'vvc'c(^1^2A3A)(£vc'c - £ee'v) 

= ^AAs <5mm3 a( Av Ac' Ai ) [a( Ac Ac' A2 )Tvc'ccc'v ( A2 Ai A) 

^ _ _ (Allb) 

+ 1 EA2A{-^ecvc'(A A2AA2Ai)rc'cvcc'v(A2AiA) 

+ (_l)Ae-A. ^_ 4,^^^,^(AiA2AiA2AA)(re'c.vee'(^2AiA) 

-a(AcAvA2)Tc 

'vcvee' 

(A2A1A))}]. 

^fv'v"w'c(^1^2A)(evv'c - evv'v") 

~ ~ (A12a) 

= 5<5aA3<5mM3[T'w'v"vv'c(AiA2A) + rvv'v"vv'c(AlA2A) 

+a(AcAv'A2) Eaja Jvv'vc(A A2AA2Ai)fvv'v''cvv'(AiA2A)], 

0^^,;„--,^(AiA2A)(£vv'c - evv'v") 

(A12b) 

= <^AA3'^MM3a(Av'Av"A2)[ia(AcAv'Ai)rvv'cw'v"(A2AiA) 
-Tvcv'vv'v" (A2 Ai A)] . 
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^w4v'c(^1^2A3A)(evv'c - Svv'e) 

= SAAsSMMsi^iTyv'evv'ci^lM^) + rw'evv'cCA-l A2A)} 
-a(AeAv' Ai){T'vev'vv'c(AlA2A) + Tvev'vv'cCAlAaA)} 

+a(AcAv'A2) Ea2a Jvv'vc(A A2AA2A1) 

x{ifvv'ecvv'(AlA2A) - ia(AvAv'A2)T;v'ccv'v(AlA2A) 
-a(AeAv'Ai)rvev'cvv'(AlA2A)}] + 

+a(AvAv'A2)yvv'vc(AiA2AA3ri)rvev'cvv'(Airi)] 

+ 5E«[(-l)^'"^^a(^vAv'Ai)/vev'vv'c(AlA2A3r2Ati) 
xT''v'evcvv'(^''"2) + a.(Av'Av'A2) 

(AiA2A3T2Au)T' 

evv'v'cv 



(A13a) 



^vv'evv'c(^1^2A3A)(£^v'c - £w'e) 

= '^AA3'^AfM3[T^cv'vev'(A2AlA) — ia(AcAv'A2) 

(A2 Ai A) + Ex,A,A a(AvAv' A2) (Al3b) 
X {a(AvAv' Ai)rcvv'evv' (A2 Ai A) 
xj;^,g-^-,(AiA2AiA2AA) - ia(AcAvAi) 
xTv 

'vcevv' 

(A2AiA)/;^,^^,,^(AiA2AiA2AA)}]. 
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^ee'vvv'c(^1^2A3A)(£vv'c " ^ee'v) 

= <5AA3'^MM3[H^ee'vvv'c(''^l-'^2A) + ree'vvv'c(A.l-'^2A)} 
-a(Ae'AvAi){reve'vv'c(AlA2A) + reve'vv'c(AlA2A)} 
+a(AcAv'A2) X^AaA J'ev'vc(A A2AA2Ai){ifee'vcvv'(AlA2A) 

-a(Ae'AvAi)(i 

evc'cvv' 

(A1A2A) + reve'cvv'(AlA2A))}] 
+<^AiA2<^A3ri5AoTee'vvv'c(Airi) + a(AeAvAi)(-l)^2 (Al4a) 
xyev.vc(AiA2AA3ri)ree'vcvv'(Airi) + ^J^a{\e' XyU) 
X {iee'vvv'c(Al A2 A3ri A'u)rvee'cvv' (uTl ) 
+/ee'vvv'c(AlA2A3r2Au)r4e,cvv'('"'r2)} + a(Ac'AvAi) 

x{a(AvAv'u)a(AcAv/A2)/cvc'vcv'(AiA2A3T2An) 

xT'e'evv'vc{uT2) + |(-l)^«-^^/eve'vv'c(AlA2A3T2Au) 
xT'e'vecvv'{ur2)}], 



^ee'vvv'c(^1^2A3A)(£vv'c - £ee'v) 

= — 5aa3 SmMs [a( Ac' Av Al ) {T^cv'cc' v ( A2 Al A) 

-ia(AeAv'A2)Tvv'cee'v(A2AiA)} (Al4b) 

+ Ew"(^eAe'Ai)a(AeA,A2){(-l)^^ 

X I'^^,^^^, (Al A2A1 A2AA)rcvv'vee' (A2 Ai A) 

+i4'v..'c(AiA2AiA2AA)r, 

'vcvee' 

(A2A1A)}]. 

388 Four-body part. 



^vv'v"v"'vv'cc'(^1^2A3A4A)(£vv/cc' - £yv'v"v"') 

(A15) 

= i(-l)^3o(A^AvAi)Fe.cvv'(AlA2A3A4A) 
X Qvv'v"v"'vcc'v' (A1A3). 
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(2) 

J^g^,^„;^^,^,j.,(AlA2A3A4A)(£vv'cc' - £evv'v") 



(A16) 



(A17) 



= 3<^AiA2'^A3A4^A0<5evv'v"vv'cc'(^1^3) + (-l)'^^a(AeAvAi) 
xFc/cvv'(AlA2A3A4A)Q^cv'v"cvc'v'(AlA3) + 
Xa(AeAvA4)Few'v" (A4A3A2AiA)(5vv'ev"cc'vv' (A4A2) 
+ 1 Swd Gevv'v''vv'cc'('"'^''^l-'^3A2A4A)Qvv'ev''vcc'v'(^^C?)- 

^ee'vv'vv'cc'(^1^2A3A4A)(£vv/cc' - See'vv') 

= 3<^AiA2<^A3A4<^Ao{Qee'w'vv'cc'(-^l-'^3) + Qw'ee'cc'vv'(-^3Al)} 
-K-l)'^H^c'cvv'(^1^2A3A4A)Qee'w'cvc'v'(AlA3) 
+ (-l)^i+^4i?;,^^,(A4A3A2AiA)Qeve'v'cc'vv'(A4A2) 

+a(AvAv'A2)a(AcAc/A4)Fcc'v'v(AiA2A3A4A) 

xQvv'ee'cvc'v'(A3Al) + a(AvAv' A3)a(AeAe' A4) 
X-^ee'v'v(A4A3A2AiA)(5eve'v'vv'cc' 

(A2A4)} +a(AeAe/Ai) 

X Su(i«(Ac'Av'<i)Qeve'v'vcv'c'('"<^) 
X Ge'ew'vv'cc' (ndAi A3 A2 A4 A) . 

389 The coefficients J, Y, Y', I, I', F, G are defined by the fohowing formulas 

^^Mp(AiA2AiA2A) = (-1)^2+^2 [Ai][A2,A2]V2 1 1 a^ \ , (A18) 

L Ai A2 A J 

>^a/3M^(AlA2AiA2A) = (-l)'^i+^^+^^+^^[Ai,A2,A]V2 

X (AMAiM, I A2M2) t 5 t } ' (^19) 

y^^^,(AiA2AiA2A) = (-l)^i+^nAi, Ai,A2]V2 

X (A1M1A1M1IA2M2) { 5 1 1 } { £ 5 t } ' (A20) 

/a/3CA^^*7(AlA2AiA2AA) = (-l)^+^i[Ai,A2,A2,A,A]V2 

X (A2M2AMIA1M1) I ^ ^ X I I ^5 1 , (A21) 

L A/x Ai J Ai A2 A J 
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TABLE IX. The expressions for phase factors Zaipip_ipi 

Aa' A/j/ A^/ Ap' Zq,'^/^/j5/(AiA2A) 

Aq A^ A^ Ap 1 

Aa A^ Ap A^ a(ApApA2) 

A^ Aq, A^ Ap a(AoA/3Ai) 

A^ Aq, Ap A^ a(AaA/3Ai) a(A^ApA2) 

A^ Ap Aq A^ a{\a\pK2) a(A^ApAi) 

A^ Ap A^ Aq, o(A^ApAi) 

Ap A^ Aq A/3 a{\a\pJ^2) 

Ap A^ A/3 Aq 1 



^a/3CAP^(AlA2AiA2AA) = (_1)A.+A,+A,+A,+Ai+Ai+A,+A+A 

X [A][Ai,A2,Ai,A2]^/2 



A Aj^ A 
A^ A2 Ai \ 
\a A2 Ai Xa 



(A22) 



F„/3Cp(AiA2AiA2A) = (-1)^^[A2, A2]i/2 { ] { t } ' (^23) 



G'Q/3Cp/iPi7a(AiA2AiA2AiA2A) = (-l)^a(AQA/3Ai)a(ApAsA2) 

= = , ,„ I -^/s Aji Ai Ai Aq, Aj7 

X [Ai,A2,Ai,A2,Ai,A2]^/^ < Ai _ A^ A2 

A^ A^ A2 A2 Ap Act 

390 In Eq. ()A22p . the last term with the brackets [. . .] denotes 12j-symbol of the second kind (see 

391 Ref.^^ (Sec. 4-19, Eq. (19.3), p. 102)). In Eqs. (|A3l)-(|A8l), the expressions for the coeffi- 

392 cients Z^' i3>p,>ii>{^i^2^)-, where {a', /3', //', P'} denote somehow permuted orbitals {a,l3^fi,v} of 

/ox 

393 the coefficient r2^^-p(AiA2A), are displayed in Tab. IIXI Particularly, we use the abbreviation 

394 ^Q/^/^/p/(AiAiO) = 2q//3/^/p/(Ai). 
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